Tracking f(R) cosmology 
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Metric f(R) gravity theories are conformally equivalent to models of quintessence in which matter 
is coupled to dark energy. We derive a condition for stable tracker solution for metric f(R) gravity in 
the Einstein frame. We find that tracker solutions with —0.361 < U) v < 1 exist if < Y < 0.217 and 
4z \nf'(R) > 0, where Y = *f% is dimensionless function, lo v is the equation of state parameter 

of the scalar field and R refers to Jordan frame's curvature scalar. Also, we show that there exists 
f(R) gravity models which have tracking behavior in the Einstein frame and so the curvature of 
space time is decreasing with time while they lead to the solutions in the Jordan frame that the 
curvature of space time can be increasing with time. 



I. INTRODUCTION 

Fourth order f(R) gravity theories can be considered 
as a candidate for solving the major challenge of cosmol- 
ogy i.e. the late time accelerated expansion of the uni- 
verse [l|, without introducing any exotic matter sources 
0(and references therein). These theories have particu- 
lar features among the other modified theories of gravity. 
f(R) modifications to GR appear in the low-energy effec- 
tive actions of quantum gravity and the quantization of 
fields in curved spacetime. These theories suggest a com- 
pletely geometric origin for both the early time inflation 
and the late time cosmic acceleration. 

These theories are conformally related to GR with a 
self- interacting scalar field Q. Although these models 
can rise to a natural acceleration mechanism, there ex- 
ists some features in them which make their viability du- 
bitable. For examples these models predict an amount 
i for PPN parameter 7, which is a gross violation of 
the experimental bound I7 — 1| < 2.3 1CT 5 Q. Albeit 
chameleon f(R) gravity models can pass the solar sys- 
tem tests but in the sense of cosmological considerations 
these theories are observationally indistinguishable from 
a cosmological constant [5|. 

Any way, our purpose here is to find out a condition for 
stable tracker solutions for metric f(R) gravity models in 
the Einstein frame. Although we pass from fourth order 
gravity to scalar-tensor gravity in which equations are 
mathematically simpler, but its physical relevance is still 
controversial [6J . However, following standard procedures 
one should not conclude equivalently about the physical 
relation between the results, ft has been demonstrated 
that passing from one frame to another can alert the 
physical meanings of the results (see Q and references 
therein). For example the stability of solutions can be 
completely different in the two frames Q. However, we 
pay our attention just to the Einstein frame. 

Consider the general action of these theories in Jordan 
frame 
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dan frame. Under the conformal transformation g^ u = 
e 2aip~^ ^ where ip = ^ In /' and prime denotes deriva- 
tive with respect to R , we obtain the Einstein frame 
action 



Se = J d 4 xy/^g[ 
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Where V = (Rf - /)/12a 2 /' 2 - We see that in the Ein- 
stein frame the scalar field couples minimally to gravity 
but couples conformally to matter fields via the function 
e~ 2av . For a spatially flat FRW space-time the modi- 
fied Friedmann equations and the equation of motion of 
scalar field are given by 

H 2 = 2a 2 (p v + p m ) (3) 



H = -3a 2 [p m (l + u m ) + Pv (l + lu v )} (4) 



(p. + 3Htp + V v — a(l - 3uj m )Pm (5) 

Where p m and p m are the energy density and pressure of 
cosmic fluid in the Einstein frame. Also, p v = ^ip 2 +V((p) 
and p v = ^ip 2 — V(ip) represent the energy density and 
pressure of the scalar field. 

In Einstein frame, the scalar field and the cosmic fluid 
satisfy the conservation equations 



p v + 3H{1 + io<p)pp = a(p{p ri 



3p„ 



p m + 3H(l + LUm)Pm = —Otlp(p m - 3p m ) 



(6) 
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and the energy density of matter p m , pressure p m , cosmic 
time t and scale factor a are related to their Jordan frame 
counterparts through @ 



Pn 



m Pm 



dt = e atp dt a 
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Where a 



and all tilded quantities are in Jor- 



It is clear from equation © that the evolution of scalar 
field is not determined only by its potential energy since 
there is a coupling to matter. Furthermore, such cou- 
plings give rise to additional forces on matter particles in 
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addition to gravity. During the matter dominated phase, 
by using equation ([5]), one can introduce an effective po- 
tential as follows 



Veffif) = V{(p)+p*e~ 



(9) 



Where p* is a conserved quantity in the Einstein frame 
[l(|, which is related to p m via the relation p m — p*e~ a,p . 



II. TRACKING SOLUTIONS 



In this section we find the condition for having tracking 
solutions. We confine our attention to the case (l v > 
which is satisfactory from the astrophysical point of 
view [ll| and the generalization to the case tl v < can 
be done with similar considerations. In the uncoupled 
quintessence model this condition requires that oj v < uj m 
and lo v be nearly constant. But here we have 

(l v = 3H(uj m — oj v )il m n v + a(l — 3u) m )ipCl m . (10) 
where 

H 2 

•f 2 = (i + u tp )p v = (i + uj v )—n v , (ii) 

Using this equation to eliminate p>, equation (|10|) be- 
comes 



V i + w« 



(12) 



At tracking era fl v <C Sl m , so we can ignore the second 
term in the above equation. In the matter dominated 
era, this shows that the condition of (l v > means that 
tp > and so we select the plus sign. Also it is not 
necessary that cdp be negative and it can takes any value 
between -1 and 1 (provided that f2 v is very small such 
that f2 v < 0.1). To obtain the tracker equation we can 
express the equation of motion of scalar field into the 
following form 
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where x 



is the ratio of the kinematic to the 



potential energy for tp and prime denotes derivative with 
respect to In a. Therefore, for a tracker solution (lu v ~ 
const) the tracker condition becomes 
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By taking the time derivative of this equation for tracker 
solution, we get 



d_ f Vfi, Say/1 + uj ip n v 
Tt [ V> ~ Q 3 / 2 1 
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Therefore, during tracking, the second term is dominated 
and we have ^(-^f) < or equivalently T < 1, where 

r is a dimensionless function defined as T = V ^ P ,Y [3 
and is nearly constant at tracker period. This condition 
is completely different from the tracker condition in the 
uncoupled quintessence model (r > 1) ^2]. By taking 
the time derivative of equation (|13| and combining with 
Friedmann equations and equation p3[) itself, we obtain 
the following relation 
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1 + u v (6 + y) 2 2(l+^)6 + y 

3(w m — ujtp) 1 — Q v XQmQtp ^ 
1+uj v 6 + y 2{6 + y) 
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Where x = f£f , y_ = x - AO m Q- 1/2 , ^ - ^ and 
-. From now we take w m = because 



A = (l-3u! m )J 



we arc considering the matter dominated era where met- 
ric f(R) gravity behaves like an interacting quintessence. 
Since fl v <C £l m during tracking, a useful relation can be 
obtained for T by expanding it as a power series of f2 v 
i.e. 



(17) 



+ [E 3 + E 4 5 + S 5 f 2 -£]Q V + 0(n 3 J 2 ). 



Where 



Ei = 6u} v /yj2(l + u v ) - 3yJ(l+LJ v )/2 - 3(1 + u v ), 

S 2 = (1 - ^8(1 + LJ V ) - yj(l + ^)/8 - 



E 3 = 18w p - 9 - (28 + 27av 
w v 3(1 + uj v f/ 2 
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(18) 
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Therefore for the tracker solution (assuming T is nearly 
constant and f2„ <C f2 m ) to the first order in u^, we have 



r 



(19) 



This equation shows that for interacting quintessence cor- 
responding to the Einstein frame of f(R) gravity, track- 
ing occurs if r < h, although there exists another con- 
straints which limit this condition. 

Another constraint which limits this interval, comes 
from the stability requirement. We require that for any 
solution for which the equation of state parameter of 
dark energy is different from tracking parameter, loq, by 
a small amount such as Su>, then Suj decays with time 
and the solution joins to the tracker solution. 



3 



III. STABILITY OF THE TRACKING 
SOLUTIONS 



Now we want to check the stability of tracking solu- 
tions with constant lo v . Consider a solution which is 
perturbed form the tracker solution, lo v = loq, by an 
amount Slo. Then we expand equation (|17| to lowest or- 
der in Slo. It should be noted that Slo and £l v have the 
same order of magnitude and so one can neglect terms 
such as VL v 8"lo or £1^/ 2 8lo. However, we do not neglect 
the terms containing powers smaller than one in f2„ and 
we apply the limit Q v — > to final solution for Slo. By 
using the equations for and expanding them to the 
first order of Slo we have = Si(wo) + oiSlo, where 



0"! 



3w 



02 = 



V 2 (l + ^o) \/2(l+^o) 3/2 ^8(1+ w ) 

LOq - 1 1 



32(l+w ) 3/2 y/32(l + LJ y 



27 

°3 = —j= (1 + w ), 



(20) 



cr 4 = 27 



1 + wo 28 + 27lo 



V8(! + wo) 



Finally, by using these equations we obtain 

5"lo + i]S'lo + £Slo = 0. (21) 



where 



r, = -[S 2 M^ 1/2 + S3(wo)], 
^ w 2 - 1 



2 v / 32(l+w ) J 



(22) 



the solution of the equation (f2"Tj) is 



(23) 



In order to have a decaying Slo , should be positive even 
if the quantity under square root is negative. One can 
easily show that r/ is positive , for any loq in the interval 
-1 < loq < 1, if fl v > 0.025. When Q v < 0.025 then the 
interval of loq for which rj is positive becomes tighter. We 
select the smallest interval corresponding to fl v ~ 0. It 
is easy to show that this interval is 



0.361 < loq < 1, 



(24) 



By taking into account equation (|19[) we obtain the final 
condition for having tracker solutions 



< T < 0.217. 



(25) 



In a similar way, the conditions such as (|25[) has been 
derived for the existence of tracker solutions for K-essence 
in [l3[ and also for quintessence and k-essence in a general 
cosmological background in [141 ]. 



The condition (|25|) should not to be confused with 
the condition T > 1 which has been appeared in the 
chameleon scalar tensor theory [13] • f(R) gravity mod- 
els can be considered as a chameleon theory for which (3 
is negative and is equal to — 1 / a/6 [l5| . In order to find 
tracker condition, the sign of (3 is important. Since in the 
chameleon scalar tensor theory, this coupling constant is 
positive, the second term in the equation (|14|) is nega- 
tive. By a similar calculation one can easily verify that 
the tracking condition is T > 1. Note that the condi- 
tion tp > is again required for having increasing density 
parameter of dark energy. 



IV. PROPERTIES OF TRACKING SOLUTIONS 

In the uncoupled quintessence model p v ~ a -3 ( 1+ "'<°) 
and always p m < and p v < 0. But here the conser- 
vation equation of the dark energy density, equation ^ , 
is not integrable even if equation of state parameter is 
nearly constant. So the tracking behavior is not clear. 
Furthermore, it is clear from ([6]) that if lo v is nearly con- 
stant then f>cp can be positive or negative. Also it is not 
obvious that whether \p v \ < \p m \ during tracking, like 
that of uncoupled quintessence model, or not. 

Here, we want to clarify these ambiguities and explore 
the tracking period with some more details. To do this, 
we can write equation ([6]) as follows 



p^ = ap v (p 



'2(1+ hy) 



(26) 



Therefore, at the beginning of the matter dominated era 
the first term is dominated and so dark energy density is 
increasing. On the other hand, when £l v becomes larger 
then the second term is dominated and so p v is decreas- 
ing. Thus, there exist a maximum at the time evolution 
of p v . We require that the dark energy density has not 
significant role in the matter dominated era and its role 
is important for us at late times. So, it is necessary to 
show that at this maximum, p v is smaller enough than 
the matter density. Let us assume that the dark energy 
density takes its maximum at the time t* , then by using 
equation (|26p we obtain 



Sir 



-13 + v//3 2 + 4(3 



(27) 



where (3 = 18(1 + lo v ). Now, by using the condition 
we have 

0.925 < n m U < 0.974. 



(28) 



Thus, the dark energy density at t* is small compared to 
the matter density. 

Moreover, it can be shown that p v decreases at a slower 
rate than p m , so t* is smaller than the time at which the 
matter and dark energy densities are equal. Firstly, let 
us to show that in this era \p v \ < \p m \- 
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After t*, by using equations (O and (JT]), the condition 
\f>tp\ < \Pm\ can be written as 



(1 



K /2 < !• 



(29) 



During tracking the density parameter of dark energy 
is small and it can be shown that the condition (f2"9")l is 
satisfied for any — 1 < u) v < 1. After tracking, in order 
to have \p v \ < \p m \, as in the uncoupled quintessence, it 
is necessary that lo v be decreasing with time. By using 
equations © and (fH?|) , it can be shown that oJ v < if 



1 

V 



dV f 



eff 



dip 



< 3a i 



/2(l+CJy) 



(30) 



after tracking Q, v is not negligible compared to 1 so 
\/2(l + uj v )/£L v ~ 0(1). Consequently, equation (|30|) 



can be written as 



r 



dV e 



eff 



d(p 



< 1. 



in which M p ; = (\/6a;) _1 . 

Now, we discuss the convergence to the tracker solu- 
tion for which lo v is nearly constant, in more detail. It is 
required that the evolution of the tracking dark energy 
density be insensitive to the initial conditions and any 
perturbation from tracking dark energy density should 
be decreasing with time. This is the main goal of intro- 
ducing the tracking solutions [12]. 

Suppose that ip is perturbed from the tracking solu- 
tion by an small amount Sip. Any perturbation in the 
"position" and "velocity" of the scalar field ip produces a 
perturbation, Sp v , in the dark energy density p v . In or- 
der to show that the fractional perturbation decays 
with time, we write the continuity equation of the dark 
energy density ([6]) as follows 

p' v + 3(1 + u v )p v = y ( 1+ 2 ^) ^/2 pm; (32) 
so a small perturbation 5p v will satisfy 



(1 + U<p) 



(n-^ - nV") s P , 



Now by using the above equations, one can write a dif- 
ferential equation for y = as follows 

Pip 



X. 



( 1 + a v) fo-V2_nVJ 



W = 0. (34) 



In order to find out fl v during the tracking era, let us 
write equation (|10|) in the form 



o. 



(35) 



0.0003 
0.0002 
0.0001 




FIG. 1: The cosmic scale factor at matter-radiation equality is 
normalized to 1. Assume that x(D ~ 10~ 3 andSl v (l) ~ 10 -3 . 



Taking into account that the equation of state parameter 
of dark energy is nearly constant for the tracker solutions 
and also f2 v <C Sim in the tracking era, this equation can 
(31) be solved analytically when the term containing ti 2 is ig- 



nored. For the sake of simplicity in solving equation (j34j) . 

3 /2 

here we will ignore also the term containing Vt^ . With 
this approximation, the above equation has the following 
solution 



6uj u 



7 a 



(36) 



Where 70 is an integration constant. By substituting this 
solution in the equation (|34[) . the solution for x( a ) is 



xw = — 
p v 



T0% /2(l + ^ y ) 



(37) 



in which xo is an integration constant and y = a 3 "W 2 . 
Also the corresponding solutions for w v — are 



In a 

x{a) = S ±L = _» a » e('- 1/4 ) 2 



(38) 



(33) This fractional perturbation has been plotted for vari- 
ous Uu,, in FigJT] We see that is a decreasing func- 
tion with time. As mentioned before, in the uncoupled 
quintessence model is constant during the tracking 
period, but here this ratio quickly reaches to zero. There- 
fore, the evolution of the dark energy density is insensi- 
tive to the initial conditions for the dark energy density. 



V. DISCUSSION 

In this paper we have derived the conditions for fourth 
order f(R) gravity models to have tracker solutions in 
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the Einstein frame. The tracker solutions exist if u v is 
nearly constant and also i) < T < 0.217 and ii) ip — 
f t In /'(E) > 0. 

For exploring the condition ii) further, let us con- 
sider only positive-definite forms of f'(R), because the 
conformal transformation — e 2av g flu is singular for 
f'(R) = 0. This forms of f(R) also are necessary in or- 
der to have positive effective gravitational coupling in the 
Jordan frame [HI]. Also, suppose that f'JR) > which 
is required for Ricci scalar stability [17] . Ricci scalar 
instability may appear if the matter-energy density (or 
equivalently the scalar curvature) is large enough com- 
pared with the matter density of the universe, for more 
details see [lU. By these assumptions, the condition ii) 
reduces to 



dR n 
- r >0 
dt 



(39) 



Therefore, although the scalar curvature in the Ein- 
stein frame is decreasing and f(R) gravity model can 
show tracking behavior, the scalar curvature of the Jor- 
dan frame can be increasing. 

Now, as an example, we want to find out a class of 
f(R) models which have constant T and so can show 
tracker behavior in the Einstein frame. Only for two 
classes of potentials the function T is constant, power 
law and exponential potentials. Exponential potentials 
lead to r = 1, which cannot pass the tracking condition. 
Thus, we consider the power law potentials. In this case, 
V(tf) takes the form 



v(ip) = v oV — = V^ n 



(40) 



Where Vo is a positive integration constant. Taking into 
account the condition i) then 1 < n < 1.28. Now, by 
using V = (Rf - f)/l2a 2 f' 2 and equation flU we get 



Rf'(R)-f(R) + U a f(R) 2 [\nf(R)Y 







(41) 



Where Uq = — 12a 2 Vo / (2a) n . Let us to solve this equa- 
tion for the case a(p -C 1, or equivalently In/' -C 1. 
By supposing that f(R) = R + ^f(R) and neglecting the 
terms containing ^>(R) with powers larger than one (note 
that n is of order unity), the equation (l4Tj) reduces to 



m'(R) - V(R) + U y'(R) n = 



(42) 



The solution of this differential equation is fi 2 (m 
and so 




FIG. 2: Effective potential of the model (J43J) , V e ff(tp) = 
Vo<p n + p*e~ av , for several values of p* . Units and constants 
have been suppressed i.e. Vo — n = a — l(note that for n = 1 
the tracker condition is satisfied). 



Taking into account the condition i) one gets m > 4.60. 
Now we show that the condition ii) is satisfied for this 
model too. When ip > then we expect that the slop of 
the effective potential which the scalar field experiences 
it is negative i.e. -j-V e ff < 0. The effective potential 
corresponding to the model (143(1 is illustrated in FigfVl 
By using equation ([9]), the condition j^V e ff < reduces 



to 



V v < ap m < ap m (t e ) 



(45) 



In this model, there exists a free parameter /i which 
can be sufficiently small to pass this condition and so, all 
condition for having tracking solutions are satisfied. It is 
important to note that the model (|43[) is similar to the 
chameleon f(R) model which has been considered before 
in the literature ||. But, in the case of the chameleon 
f(R) gravity models, theories of the kind (H3")) are com- 
patible with observation in the range of the parameter 
< m < 0.25 [l5j . Therefore, this class of chameleon 
f(R) gravity models cannot pass the sufficient condition 
for having tracker solutions. 
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f{R)=R-H 2 (\-m) 



(43) 



Where m = — ^ and 

n—l 



-Uo 



[n m (l - m) 2 }~ 



(44) 
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